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Abstract 

We calculate the one-loop correction to the soft-gluon current with massive fermions. This current is process 
independent and controls the singular behavior of one-loop massive QCD amplitudes in the limit when 
one external gluon becomes soft. The result derived in this work is the last missing process-independent 
ingredient needed for numerical evaluation of observables with massive fermions at hadron colliders at the 
next-to-next-to-leading order. 



(~| . 1. Introduction 

The main obstacle for the numerical evaluation of collider observables at higher perturbative orders is 
the presence of infrared (IR) (soft and coUinear) divergences in parton level calculations. These divergences 
r^ I cancel in observables, but need to be regularized in all intermediate calculations by introducing an appro- 

priate parameter (typically dimensionally) . It is the need to keep track of such a rcgularization parameter 
that prevents the ab initio application of straightforward methods for numerical integration. 
^ i At the next-to- leading order (NLO), this complication can be evaded within the so-called subtraction 

^^ ' method. Its basic idea is simple: first, one utilizes the universality and factorization property of IR singu- 

OO . larities to construct an approximation to the corresponding real emission partonic amplitude. This approx- 

imation is simple enough to allow the analytic extraction, and eventually cancellation, of IR singularities. 
Second, one explores the fact that the difference of the full amplitude and its approximation is IR finite 
r^ \ and therefore can be integrated numerically in a straightforward way. This approach, in effect, splits the 

task of performing a complicated divergent integration in two: first, a divergent integration of a simpler 
quantity and, second, a complicated, but finite numerical integration. The subtraction approach can be 
applied to processes with massless and massive fermions. Several subtraction schemes have been proposed 
[l|, I2, y, y, la, |6[ and have been successfully used in a large number of applications. 
^^ . To construct a subtraction scheme at next-to- next-to-leading order (NNLO) one needs to know, among 

5— ( ' others, the limiting behavior of one-loop amplitudes when one of the external on-shell partons - a gluon - 

^ becomes soft. For the case of massless fermions (like massless QED or QCD) this problem has been studied 

in specific cases in Refs. [7|, [Sj and later in Ref. |9]. The goal of the present work is to generalize the process- 
independent approach of Catani and Grazzini 9] to the case of massive fermions. With the result of one of 
the authors for the treatment of double real radiation jlOi, TlJ, the result derived in this paper represents 
the last missing ingredient for the construction of NNLO observables with massive fermions. Applications 
are top and bottom (and charm) production at hadron colliders, deep inelastic scattering and processes at 
lepton colliders. 

This paper is organized as follows: In Section [2] we present the factorization of one-loop amplitudes in 
the soft limit and introduce the one- loop soft-gluon current. There we also present the derivation of the 
one-loop soft-gluon current in terms of a set of scalar integrals. In Sections [2.1112.21 and 12.31 we present the 
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explicit result for the UV unrenormalized one-loop soft-gluon current for all phenomenologically relevant 
kinematical configurations. In Section 12.41 we discuss the properties of the results and explain a number 
of checks we have performed. In Section |2] we derive the soft limit of squared matrix elements as needed 
in specific applications. In Section |4] we describe the UV renormalization of the bare one-loop soft-gluon 
current followed by a summary. We have added a number of appendices containing most of the technical 
details. The evaluation of all scalar integrals is detailed in [Appendix A[ In [Appendix B| and [Appendix C] 
we independently derive two limiting results for the UV renormalized one-loop soft current: its small- mass 
limit and its pole terms, respectively. Finally, in [Appendix D[ we discuss the analytical continuation of the 
bare one-loop soft-gluon current evaluated in different kinematical configurations. 

2. Amplitude factorization in the soft limit 

Consider the amplitude Ma{n + 1; g) for producing n+ 1 on-shell partons. [^ Let at least one final state 
parton be a gluon, and let a = 1, ... , N^ be its color index and q its momentum q^ = 0. It is useful to think 
of Ma{n +\]q) as a wide-angle scattering amplitude, i.e. all kinematical invariants formed from its external 
momenta are large. The structure of such amplitudes is very well understood through at least two-loops 
IS [13, [3 H lli, [13, [IS H 0, il ii, [11, [11 [IS S [13 in both the massive and the massless cases. We 



are next interested in the limit when the external gluon becomes soft q — ^ 0, or more precisely its momentum 
scales as: 

g-> Aq, A^O. (1) 

Clearly, in the limit ([IJ, Ma{n + 1) is not a wide-angle scattering amplitude anymore. Instead, it satisfies 
the following factorization property: 

Ma{n + l;q) = Jaiq)M{n) + 0{X). (2) 

In the above equation, the amplitude M{n) is the wide-angle scattering amplitude obtained from Ma{n+\] q) 
by removing the external gluon with momentum q, and Jq {q) is the process-independent soft-gluon (eikonal) 
current whose derivation is the main goal of this work. We have adopted a concise notation in Eq. ^, but 
have made explicit the dependence on g, a and n. 

Each one of the factors in Eq. ^ depends on the dimensional regularization parameter e = (4 — d)/2 
and has a loop expansion in powers of the strong coupling constant through one loop: 

■Uq) - gst^'{j^^\q) + Jl^Hq) + .-) , 

M{n) = M^°'> (n) + M^^'> (n) + . . . , 
Ma{n + l;q) = Aff )(n + 1; q) + Af «(n + 1; g) + . . . , (3) 

where the dots stand for terms at higher orders in as- The notation adopted in Eq. ^ does not make 
explicit the powers of the strong coupling as- For example, Ja denotes the leading order result for the soft 
current (given explicitly in Eq. Q below), Ja stands for its next-to- leading order in as, and so on. The 
reason for choosing this notation is that the leading order amplitude M'-^^n) contains a process dependent 
power of the strong coupling constant. Thus, our notation reflects the only relevant information: the power 
of the strong coupling relative to the leading order amplitude M^'^^n)- 

Our considerations apply for both bare and UV renormalized amplitudes. For now we consider bare 
amplitudes and will return to the UV renormalization in Section [51 

Both in the massive {pj > 0) and massless {pf = 0) cases the tree level soft-gluon current reads: 

n fj. n 

Jr\ci) ^Y.Tt^^^Y.Ttet . (4) 

z=l P^ ■ 9 i=l 



^Unless we state otherwise, we do not make a distinction between initial and final state partons. 
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where Ja (q) = e^{q)Ja (q)- Throughout we follow the conventfons of Ref. ^] for the signs of color 
generators. The one-loop UV un-renormalized soft-gluon current Ja (q) reads: 

n 

J^^'\q) = ^fabc Y. T'T^ i< - <) 9S\e,q,p,,p,) . (5) 

For the calculation of the one-loop soft-gluon current in Eq. ([S|) we follow the strategy developed in 
Ref. [9[. The approach consists of the evaluation of all one- loop diagrams connecting (on-shell) external 
legs and attaching a real gluon to either the external legs or the gluon propagator (the virtual gluon). The 
calculation is performed in the cikonal approximation, i.e. the real and virtual gluons are treated as being 
of similar magnitude, and energy-momentum conservation is enforced. 

Following the terminology introduced in Ref. [9|, we split the results in IP and 2P contributions. The 
IP contributions arc defined as the ones that depend on a single external hard momentum pi, as opposed 
to the 2P contributions that involve two hard momenta pi and pj. In the following we calculate the 2P 
contributions and show that they are separately conserved; then, adapting the arguments given in Ref. [9|| 
one can show that the IP terms do not contribute to the soft current. 

Our starting point for the calculation of the 2P contribution JJ^Lp) to the one-loop soft-gluon current is 
the sum of the three diagrams 4(a, 6,c) given in Ref. [91]. g We neglect all scaleless integrals. The sum of 
diagrams is gauge invariant as also explained in [9|. The term ~ k ■ e{q) is reduced to scalar integrals. The 
reduction differs from the one in Ref. [9| since in the case at hand, it produces terms that explicitly depend 
on the masses ml :. Applying partial fractioning and omitting scaleless integrals, we arrive at the following 

expression for the function g^ ■ in Eq. ^ : 

' ml (j>, ■ g)2 - 2(p, • p, ) {p, ■ q) {pj ■ q) + m] [p, ■ q^ 

X \{Pi-(l){P]-(l) [P] ■ q)Mi + [p, ■ q)Mi 

+ 2 (Pi ■ l) [(P' ■ Pj)iP^ ■ l) - ^UP] ■ q)] M2 + 2^Pt- <l) [(k • Pj)(Pj ■ Q) - ™i(K • <?)] M2 

+ [ip, ■p,)ip, ■ q){p, . q) - mKp, ■ qf - rn^P. • qf] ^^^^^^^^^1^ M,] . (6) 

Pi ■ Pj 



The bare coupling a'g = a^S'e/(27r) with Se = (AttY eyip{-ejE), Mk = Mk{pi ^ Pj), k = 1,2,3 and the 
integrals Mi^2,3 can be found in [Appendix A[ Noting that M3 = M3, it is apparent that gj - — 9ji ■ From 
Eqs. (I4I5P it is evident that the massive soft-gluon current is conserved through one loop. This follows from 
color conservation (as explained in Ref. j:9j) and the identity q ■ ei = 1. 

Next we present our main result, namely, the explicit expression for the function g\j ■ There are three 
kinematical regions for which this function needs to be computed. In all three cases we consider pf = ml > 
and take pi, as well as the momentum q of the soft-gluon, to be in the final state. Thus, the three kinematical 
configurations are defined as: 

1. p^ = 0, Pj incoming, 

2. p^ = 0, Pj outgoing. 



3. p^ = m^ > 0, Pj outgoing. 



"3 

d2 _ „„ 
3 3 

We do not have in mind phenomenological applications with massive quarks in the initial state, but for 
completeness, have calculated and presented below all required ingredients for such applications as well. 



^Since, for these diagrams, at the integrand level the eikonal approximation is identical in the massive and the massless 
cases, we can simply use the sum of the expressions given in Eq. (46,47) of Ref. Q . We have verified the agreement. 



2.1. Case 1 

We have evaluated this kinematical configuration directly, as described above, by substituting the explicit 
results Eqs. (|A.2[|A.4llA.14|) for the scalar integrals into Eq. ^ and then expanding in epsilon to the desired 
depth. The result for the un-renormalized one-loop soft current reads: 



g^^\Case 1) ^ R^'^ + zttI^^ 



j[Cl] 



[Cl] 



b I 2(p, ■pj)n^ 



2{p, ■ q)2{p, ■ q) 



E^"( 



R(f\Cn_^^^jin)lCl 



The bare coupling a^ is introduced in Eq. ^ and: 



r(-2)[Cl] 
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(-i)[ci] _ 



1 

2mf{pj ■q)\n 



= 4 [{p^ ■ pj){pi ■ q) ~ mjipj ■ q)] U2 (l " y) + ^liVj ' <?) In^ (y) 
2 -2(P» • P]){Vi ■ q) + m1{pj ■ q) 
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(-2)[C1] 
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(-1)[C1] 
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Rs R. 



(i)[ci] 



Rs R 



(2)[C1] 



Li.s 1 - 



Lis 



-0 



= 4 [{p, ■ pj ) (pi ■q)-m^ {pj ■ q)] 

+2 [ip, ■ p,){p, ■ q) - mi{p, ■ q)] In (l - ^) In^ (^) + ^ 



40 (k • Pj)(p» ■ q) - 26m|(pj ■ q) 



™f(Pj -9) j^3 /^a 



(t) 



1 
2' 



-4(p» ■ P]){pi ■ q) + "T-f (pj • q) 



+ 4 [(Pj • Pj){Pi ■ q) - m^iPj ■ q)] Li2 (1 - 



= 0, 

= m^(pj-q)ln^ 



2 5(2(Pi • P])iP^ ■ q) - mjiPj ■ q)) 



4 [{Pt ■ Pj){Pi ■ q) - mfipj ■ q)] Lig (^y j - C 

+2 [(p, . p,){p, ■ q) - m,2(p, . q)] In (l - ^) In^ (^) 



4 [7(pi • Pj){Pr ■ q) - 5m^{pj ■ q)] 



-In 



(t) ' 



'2{pi ■Pj){pi ■ q) - 5m1{pj ■ q) 



+ 4 [{Pi ■ Pj){Pi ■ q) - rnl{pj ■ q)] Lij (l " y) j . 



'4 [{Pi ■ Pj){Pi ■ q) - ml{pj ■ q)] 



., ,_±Uo.(:-|)-L,(|) 



4 458(p, • pj){p^ ■ q) - 2Umf{pj ■ q) 



-TT 

In ('^1(4 



720 



[{Pi ■ Pj){Pi ■ q) - mfipj • q)] Lig (^1 - y j - 2C3 [2{pi ■ Pj){pi ■ q) - m^,{pj ■ q)] j 

in (1- I) ...'(f) 



; -4(P» • P]){Pt ■ q) - mUPj ■ q) ,2 /" 



(I) 



,{Pt-Pj){Pi -q) -mf{pj -q) 



-2{p, ■ pj){pi ■ q) + 3 m1{pj ■ q) ^^4 (a 



12 



(y) " '^^y [(P« ■ Pj)(P' • 9) - m^iPj ■ q)] Li2 (l - y 



We find it convenient to express the result through the variables Rs and ai defined as: 

m1{pj-q) _ m]{pi-q) 



Rs =4: [mj{pj ■ q) - 2{pi ■ Pj){pi ■ q)] , ai = 



{p^ ■q){pi -Pj) 



, Uj 



{pj ■q){pt -Pj)' 



(9) 



We introduce the variable aj for later use. The result in Eq. ([5]) is also available in electronic form. 



2.2. Case 2 

To derive this case we have performed an analytical continuation of the result from Case 1. When p^ = 
the continuation amounts to exchanging pj — >■ —pj (see [Appendix D[ ) . It is easy to see that Eq. ([5]) remains 
unchanged under this transformation, i.e. the result in Case 2 is identical to that in Case 1. 

2.3. Case 3 

We calculate the result for Case 3 up to and including terms of 0{e). We note that the 0{e'^) term 
contributes only, if multiplied by a term ~ 1/e^ originating from the phase-space integration. Such a leading 
pole may only be due to the emission of soft and coUinear radiation. Both partons i and j beeing massive, 
collinear singularities are regularized and do not lead to poles in e. 

We compute the result for the one-loop soft-gluon current in this kinematical configuration directly, as 
expansion in e. Details of the calculation can be found in [Appendix A For completeness, and with more 
formal applications in mind, we have also calculated all integrals in the spacelike region, where the massive 
momentum pj is incoming; see Appendix A[ The relation between the results in the two kinematical 



configurations is discussed in [Appendix D 



The explicit result for the soft-gluon current in the kinematics of Case 3 reads: 

r[C3] b I 2(p, ■pj)fM^ " ' 
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Qs 
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V \lb \aj J 4 ^ ' 
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2 
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1 - 



47r" 



-4 (m^^(p, • qf - mlipj ■ qf) ^ (2 ln(w) + In {x^ + l) + ln(a;)) In ( — 

+ {m]{Pt ■ qf + mlipj ■ qf - [p, ■ Pj){p, ■ q){pj ■ q)) [^ 

321n(2) (-ln(a, +U + 1) -\n{aj + v + 1) - 2 In (x^ + l) - ln(x)) +641n2(2) 
-hl6 In (a;^ -|- l) (ln(a., + v + I) + h\{aj + w + 1)) + 16 ln(a., + v + I) h\{aj + v + I) 
+161n(x)(ln(-aj + u + 1) + ln(aj + w - 1) + 21n (x^ + l)) - 161n^(x) 



(10) 



(11) 



+81n ( — ) (ln{aj + v - I) - H-a^ + w + 1)) + 41nM — 
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2a, 
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2u + 2 
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16Li2 



2a, 



V + aj + 1 
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(-1)[C3] 



1 
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2w ^ ^ V24 12i; 



- ln(u) ln(x) + - ( 1 + - ] ln(x) ln(l + x^^ 
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i'^'^'-^'^-^s 



(m|fe.g)2 + m,^fe.,)2)ln2(^^) 



R 



(1)[C3] 



+4 (mj2(p, ■ g)2 + m2(pj ■ qf) W{x) - 4 

i (- ln(«) (ln(x) In (x^ + l) + n') + i^ + ^ 
w \ 12 2 



m^^iPt-q)^ -m1{pj -qf 



In ( — ^ 1 ln(x) 
a. 



'"wU'-'d 



Lis (a;2) _ 3 
4 



- - In^ (a;2 + l) - ^^ 

A V / 0/1 



5^ 
24 



^^ + ^ ) In (x^ + l)-\ (2Li3 (1 - x^) + Li3 (x^)) 



1 r 



^/321n^(x) 280C(3) „„, , n ,x, ,/ x 
fe • P,)(P» • q){Pi ■ q) ( 7^ ^ - 32 In (x^ + l) ln2(x) 



2^"M\ , /..,„2,.2 , ,^ 407r2^,^^^^^ 



3 3 

+ ln (x2 + 1) f 367r2 - 41n2 ( ^)) + (48 In^ {x^ + l) - 
-fln3(x^ + l)) 

(m2(p, . g)2 + m2(p, . g)2)( (31n (x^ + l) + ln(x)) In^ (^^ 

70 
+281n^(x2 + l) -441n^(x2 + l)ln(x) - —n^ In {x^ + l) + 281n (x^ + l) In^(x) 

o 

28, 3, , 2 2w X 224C(3) 
-— ln3(x) + -^2 ln(a;) + ^ 



{m](p, ■ qf - ml{pj ■ qf) ^ ( a_ 

a 



^ "' ' In ( — ) ( 4Li2 (x^) + 4 In (x^ + l) ln(x) + 8 ln(u) ln(x) 



+ (™?(Pj -qf +m]{p, -qf - {p, ■ pj){p, ■ q){pj ■ q))(^ln^ i-^j 



-InM — I21n(«) 



In^ 



2(ln(Q;j +v + l) + Ini-aj + v + I) + ln{aj + v - 1) - 31n(2)) 



147r2 



In ( — ) (21n2(w) - 12\n^{a^ + v + l) + 6\n^{aj + w + 1) - 41n2(x) -121n2 C^^ + l) 
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+28M2)ln{a, + V + I) + Aln{v)iln{a^ + V + I) -21n(aj + -y + 1) + ln(2)) - 10111^(2) 

-41n(ai +V + 1) ln{aj + v + 1) - 81n(2) \n{aj + v + 1) 

+8(2 \n{ai + v + 1) - \n{-aj + v + 1) + ln{aj + v - 1) + \n{aj + v + 1) - 3 ln(2)) ln(x) 

+24(ln(2) - ln(aj + v + 1)) In (x^ + l) + 24 ln(a;) In (x^ + l) - r^r^) 

32 

+— \n^{a, +v + l) + 121n^(aj + w + 1) - 32 ln(2) ln^(aj + v + 1) - 81n(w) ln^(aj + -y + 1) 

4 
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+56 ln(2) Ma^ + v + 1) ln(.T) - 8 \n(a^ + v + 1) ln(aj + v + 1) ln(a;) - 36 ln^(2) ln(a;) 
+16 ln(2) ln{aj + v + 1) ln(a;) + 32 In (a, + w + 1) In (.t^ + l) + 80 In^ (2) In {x^ + l) 
+32 In^iaj + w + 1) In (a;^ + l) - 80 ln(2) Ma^ + w + 1) In (a;^ + l) 
+161n(a, +V + 1) ln{aj + -y + 1) In (.x^ + l) - 801n(2) ln(aj + w + 1) In (a;^ + l) 

+ 16(-4 ln(ai + v + 1) - ln{a^ + -y + 1) + 5 ln(2)) ln(a;) In (.t^ + l) ^ ' ' 



+41n(a, +u + l)ln^(Q:j +w + l)-401n(2)ln^(aj + v + 1) - 81n(i;) In^(x) - -tt^ ln(i;) 
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The polynomial Qs, the "confornial" variable x and relative velocity of the quark pair v read: 
Qs = 16 (m^jiPt ■ qf ~ 2{p, ■ Pj){pi ■ q){pj ■ q) + mfip-j ■ qf) , 



X = V(l-^')/(l + t'), 



9 9 

m~m- 



'1--^^. (12) 

Beyond order 0(6"), the result for the one- loop soft current cannot be expressed in terms of standard 
poly logarithms. Multiple polylogarithms appear as evident from Eq. (jA.131) . At order ©(e^) we have 
combined all functions that are outside the class of the standard polylogarithms into the function: 

ri ln(l-i)ln(l-t^) 

Fcixi,X2)^ dt ^ "-^, (13) 

■Jo ^ ~* 

which can be expressed in terms of multiple polylogarithms of weight 3, see Eq. (B.22) in Ref. [28|. The 
result Eq. (|TT|) is also available in electronic form. 

2.4- Properties and checks 

The purpose of the overall d-dimensional prefactor in Eqns. (jTllOp is to extract exactly the leading power 
scaling behavior of the one- loop soft-gluon current in the limit g — >■ 0. The remainder is given as expansion 
in e, which has a well defined limit g — 5- 0. This is easy to see since it is invariant under independent rescaling 
of the momenta q,pi,pj. 

The result for the one- loop soft-gluon current satisfies a number of consistency checks. Eq. (fTU|) has a 
well defined limit, when either one of the masses rui or ruj vanishes. In the limit nij — > 0, it agrees with 
the result for the soft current in Case 2, as it should. Note that this agreement is a non-trivial check on 
the analytical continuation used to derive the result in Case 2 from that in Case 1. We have numerically 
checked the result for the hardest integral M3 in the "time-like" kinematics Case 3 (see [Appendix A\ . 



We have verified that the soft current has the correct behavior in the small mass limit (see Appendix B| 
for details). The massless limit rrii = 0, rrij = of the one- loop un-renormalized soft current is regular, and 
the results for the soft current in all kinematical regions reproduce the massless results of Ref. [9] . 

We have also verified that the pole terms of the one-loop soft current agree with what is expected based 
on the structure of the singularities of massive gauge theory amplitudes (see [Appendix G\ . 

3. Squared matrix elements 

The knowledge of the soft-gluon current makes it possible to construct an approximation to the squared 
one- loop matrix element for any process in the limit ([T]). As indicated in Eq. ([2]), this approximation is 
correct up to power suppressed terms. The result Q for the tree-level current is exact in e. The one-loop 
current ([5]) is calculated as an expansion in e which is deep enough to allow the derivation of the terms 0{e^) 
in any observable at NNLO. 

In the limit ([ij the square of a Born amplitude reads: 

(MW(n + l;g)|M(0)(n + l;g)) = 

{71 n \ 

Y, e,,(M(o)(n)|T, • T,\M("\n)) + ^ C,e,,(M(°)(n)|M(o)(n)) + 0(A) . (14) 

Above we introduced e^j = e^ • ej and Ci = Ti ■ Ti is the quadratic Casimir appropriate for the parton i. 



The interference term between the Born and one-loop amphtude in the limit ((Ij reads: 
{Mi'^\n + l;q)\A4^\n + l;q))+c.c. = -iTrasfi^'l 

n n 

2Ca Y. (^^J- - e^) R^AM^"Hn)\T^ ■ T,\M^°\n)) - 47r ^ e,kI^AM^"Hn)\r''-TtT^T^\M^"Hn)) 
J2 e,,(MW(n)|r,.r,|Af(i)(n))+c.c.| + ( f^ C,e,,(Af (")(n)|Af (i)(n)) + c.c. J I + 0(A) , (15) 

^i^j=l J Vi=l / J 

where we have split g^^ = Rij + inlij into its real and imaginary parts to be found in Eqns. (J7ll0p . 

4. UV renormalization 

Up to here we considered bare amplitudes. In practical applications one works with UV renormalized 
amplitudes. The UV renormalized one-loop soft-gluon current is very easy to derive. One needs to recognize 
that though that loop order no mass renormalization enters. Therefore, all one needs to do is coupling and 
field renormalization. Since we consistently set to zero scaleless integrals, the only correction one needs to 
take into account is self-energy contribution in the soft-gluon leg due to the massive flavors. It will be most 
convenient to work in a scheme where all massive flavors are decoupled, i.e. the coupling is running with 
Uf light flavors only. Then the heavy quark loop contributions into the external gluon leg will be canceled 
by the decoupling correction. Therefore, in order to obtain the UV renormalized current from the bare one, 
one only needs to perform coupling renormalization a'gSe — as(l — /3oas/(27re) + Olog)) in the first line of 
Eq.Q. This procedure amounts simply to adding the term ^ (3o in Eq. (IC.SP to the bare current; see also 
the discussion following Eq. (jC.Sp . 

5. Summary 

In this paper, we have studied the behavior of one-loop QCD amplitudes with an arbitrary number of 
massive fermions in the limit when one external gluon becomes soft. Similarly to the well known massless 
case, we find that in the limit ([T}, any amplitude factorizes, up to power suppressed terms, into a product 
of a simpler amplitude and a process independent function: the soft-gluon current. 

We have explicitly calculated this current through one loop. This result enters the evaluation of any 
cross-section with massive fermions at next-to-next-to-leading order within a subtraction approach. An 
immediate application for this result is the calculation of the ti cross-section at the second perturbative 
order. 

We have performed a number of non-trivial checks on our results. We have verified that they correctly 
reproduce their small-mass limit and pole terms, that we independently predict. We have performed a 
number of numerical checks on the non-trivial integrals. 

The explicit result for the one-loop soft-gluon current with massive fermions is much more complicated 
than in the massless case. While it is possible to derive an exact result valid in d-dimensions, we have 
explicitly presented the final result in a form suitable for practical applications: we have extracted the 
current's leading behavior exactly in d-dimensions, and expanded the rest in e in a form appropriate for 
calculating observables at next-to-next-to leading order. The functional form of the result is significantly 
more complicated and involves multiple (Goncharov) polylogarithms. 

As a by-product of our calculations, and as an additional crosscheck, we have worked out the analytical 
continuation from space-like to time-like kinematics for a multiscale problem. In higher orders in e this 
procedure involves multiple polylogarithms and opens an interesting subject for further investigation. 
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Appendix A. The scalar integrals 

The one-loop soft-gluon current can be expressed through the following integrals: 

d'^k 1 



Ml = $ 



i{2Try 
d'^k 


[P 


][{k 


1 


-Pr 


k] 




i{2TTY 

d'^k 


[P 


][P^ 


■ k+pi 


1 


-Pj 


■ k] 



^' " * ' ^{2^Y W] [{k + qY] [p^■k+p,■ q] [-P, • k] 

where each propagator has an implicit +i5 imaginary part. The momenta Pi^Pj can be massive or massless 
and the momentum g, corresponding to the soft-gluon, is assumed outgoing and massless. The normalization 
factor is $ = 87r2(47r)-'e''^^. 

The simplest integral to evaluate is Mi: 

Next we consider the integral M2. Its full q dependence can be extracted, and the remainder expressed 
through a one-dimensional integral: 

f 1 ~'^' 

dtr'^Ufmf + {l~tYnY^~2t{l~t){p,-pj)-iS'> . (A.3) 



X 



This one-dimensional integral can be evaluated in terms of 2-^1 -type hypergeometric functions. After some 
rearrangements and using standard relations between the hypergeometric functions we obtain: 

A/2 = $ ^^r(-6)r(26) [-{p., ■ q) - tS]-^' [-2(p, . p,) - zS]-'+' 

[ r(l + e)r(l-2e) _i+2e.-e 2/3,f / 1-v 

^ i TT^ ^ " Pj -7^2/^1 1,1- e,l + e;-— — 



We have introduced /3fe = m\/{—2{pi -pj) — iS), k — i,j, and the relative velocity v defined in Eq. ([T^. The 
hypergeometric function can be expanded in series in e to any desired depth with the help of [29]. 

The most complicated integral is M3. We first apply Schwinger a-parameterization. The two integrations, 
corresponding to the two propagators quadratic in k, can be transformed in the usual way: 

daida2 = / ada / dy , (A. 5) 

^0 Jo 

and the integration over a performed. 

In order to extract the scaling behavior of the integral in the limit g — s> 0, we rescale the a-parameters 
a3,4 corresponding to the two propagators that are linear in k: 

cta ^ a3/\pi-q\ , ai-^ ai/\pj -ql. (A.6) 
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We note that the invariants {pi ■ q) and (pi ■ pj) are non-zero, although their signs change depending on the 
kinematical configuration. Next we change the variables 0.3^4 along the lines of Eq. (jA.5|) and perform the 
integration over the infinite range, arriving at the following two-dimensional representation for M3: 



-2+£ 



M. 



$ 



16 



-r(-e)r(2 + 2e: 



} f dt f dy (t^Y^^ 

\pi-<im -il Jo Jo \ (Pi-Q) 



+ a-t) 



2 "^3 

(PJ ■ 9)^ 



-2i(l - t)-^^^p^ - i5 ) {-ya, - ta, + (a, + a,)ty - zSy^-^' 

\pt ■ q\\P] ■ <i\ 



(A.7) 



where ak = (jpk ■ q)/\Pk ■ q\ = ±1- Note that the signs atj also determine the sign of (pi ■ Pj). 

The evaluation of Eq. (jA.7|) is hardest in the phenomenologically relevant "time-like" (TL) kinematics 
of Case 3, where Ci = ctj = 1. Performing the y- integration we get: 



MP)=$ '' 



+ {l-tf 



lQ{pi -pj 

m]{pi -q) 
(pj ■q){pt -Pj) 



-r(-e)r(l 4- 2e) (-1 - i(5) 



-2t{l-t) -i6 



isy 


2e 




Pr ■ Pj 




\iPi 


q)iP3 


q) 


f-i 


-2e 


-(1 


-t)-' 


-2£ 



771 



1+e pi / 

\ {pt ■ q){pt ■ Pj) 



(pj ■ q) 



i\f3 



l-2t 



(A.8 



We evaluate the above integral as expansion in e using the formula: 



_. 00 



k=0 



e" (\r{z) 



Extracting the imaginary parts is particularly laborious since both roots x\ 2 of the polynomial, which is 
quadratic in i, are inside the integration range: < x^ < 1/2 < x\ <1. The expressions for the two roots 
in the time-like case read: 

(A.9) 



Uj +\ — V 



1 ^2 



Uj +1 + V 



riTL) ■ 



with ai , aj defined in Eq. ([9]) . The explicit result for the integral M3 is rather lengthy. It is supplied in 
electronic form. 

The integral M3 in Eq. (IA.7P is much easier to calculate in the "space-like" (SL) kinematics ai = — CTj = 1; 
in the following we present its derivation for completeness. Performing the y-integration and after some 
simplifications we obtain: 



M^^^^ = $ 



-2 + € 



-r(-e)r(i + 26) 



Pt ■ P3 



m: 



V {pi ■ q){pt ■ P3) 



\P]-q) 



Hi-tf^"^^"^ 



lQ{pi-Pj) \{pi ■ q){pj ■ q) 

+ 2t(l - t) ) [(1 - i)-l-2. ^ (^_i _ ^^y2c^-l-2c^^ 



{pj ■ q){p^ ■ Pj) 
The one-dimensional integral can be evaluated in terms of the Appell hypergeometric function Fi : 

.-2+. r(-e)r(2e) 



a4^^^ = -$ ^ 



Pt ■ P] 
16 (pi ■ q){pj ■ q) V ipi ■ q){pj ■ q) 



(A.IO) 



(A.ll) 



(-1 - ^S)-'^ a)F, [-2e, -., -., 1 - 2e; ^, ^ j + atF, (^-2e, -e, -e, 1 - 2e; ^, ^ 
In the above equation we have introduced the following notation: 



X, 



.T, = 



a 



3 



I — V aj — 1 + V 
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(A.12) 



The quantities xf 2 are the two roots of the polynomial quadratic in t appearing in Eq. (jA.10[) with x^ < 
and a;| > 1 and a^ , aj are defined in Eq. ^ . 

The Appell functions can be expanded in e in terms of multiple polylogarithms l^imk,....mi{tk, ■ ■ ■ ,ti) 
with the help of the library Nestedsums [30[ (see also |3l|): 

Fi(-2e, -6, -6, l-2e;t,y) = l + e^ [2Li2(t) + 2Li2(y)] 



4Li3(i) + 4Li3(y) - 2Si,2{t) - 2Si,2{y) - 2Lii,2 [ -, y 

y 



2Lii,2(^,t 



+e'' 



8Li4(i) + 8Li4(y) - 4^2,2(0 - 4^2,2(2/) + 25i,3(i) + 25i,3(y) - 4Lii,3 ( -, 2/ 



4Li,3(f,.) 



+2Lii4,2 f 1, -, 2/ j + 2Lii4,2 (l, |, t) + 2Lii,i,2 f -, 1, y j + 2Lii4,2 (|, 1, i) 



-2Lii,i,2 ( -^7 -J 2/ 



' J- ' 

y t 



0{e') 



(A.13) 



The functions Li„(t) are the usual polylogarithms and Sn,p{t) are the Nielsen's generalized polylogarithms. 
We follow the conventions and definitions of Ref. [30| (see also Ref. [32|)- The numerical evaluation of 
multiple polylogarithms has been automated in Ref. [33|. 

The results given above are sufficient to explicitly derive the one-loop soft-gluon current in the kinematics 
where one of the massive quarks is in the initial state. Such formal result is of interest, for example, in 
studies of the properties of massive gauge-theory amplitudes. 



(SL) 

In the case of one non-zero mass, as needed for Case 1, the integral M3 reads: 

^-2+^r(-e)r(l + 2e) f 2p,-p, 



M. 



(SL), 



|t?i^ — 



=n = $ 



16 [pi ■ q){pj ■ q) V (P» • l){P] ■ 1) 



(A.14) 



r(l + e)r(-2e) 

r(i-e) ■ 



Fi 



-e, 1 + 6,1 -e;l 



2 



{-l-i5) 



-2e 



2F1 



(-£,-£, 



l-e;l 



2 



In the massless case, the result is simple and agrees with the one given in Ref. [9[. 

Finally, we remark that even in the case of equal masses, the result for the one-loop soft-gluon current 
is a function of two independent parameters, i.e. the case of two un-equal masses is not more complicated 
than the equal mass case. This is evident, for example, from Eq. (jA.lip . 

Appendix B. Small-mass limit of J^'^^^ 

Following the methods of Ref. [15| one can independently derive the leading behavior of the one-loop 
soft-gluon current in the small-mass limit. Considering UV renormalized amplitudes, assuming all non-zero 
masses are equal and then taking the small- mass limit of both sides in Eq. ([2|), one easily derives: 



Ja{q]rn ^ 0) 



^(m|0) 



Zy'r'Ja{q;m^0) + O{m^). 



(B.l) 



In the above equation we take the strong coupling running with til + I flavors, i.e. the heavy flavor is active. 



The factor Z^"^' ' is given in the appendix of Ref. [15[ through one loop and to all orders in e. One can also 



check that upon decoupling the heavy flavor (in d-dimensions) the Z-factor in the above equation is exactly 
compensated through one loop, i.e. in the decoupling scheme, Eq. (jB.ip simplifies to: 

Ja{q\ m ^ 0) = Ja{q\ m - 0) + ©(m^) . (B.2) 

The decoupling relations in d-dimension can be found, for example, in Ref. [3J|; see also Section 21 
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Appendix C. Pole structure of J^^^^ 

One can provide an independent derivation of the poles of the UV renormalized, one-loop soft-gluon 
current from the known structure of the singularities of one-loop gauge theory amplitudes. Considering the 
amplitudes appearing in Eq. ^ as wide-angle scattering amplitudes, we can decompose them into jet, soft 
and hard functions respectively |12l |: 



Main + l;q) 
M{n) 



I X Sab ■ Hb , 

i X a ■ h . 



(C.l) 



The jet functions /, i are diagonal in color. At one-loop they contain double and single poles. The soft 
functions S, a are color matrices that have single poles only. The hard functions iJ, h are finite color vectors. 
Applying the decomposition (|C.ip to Eq. ^ and expanding each factor through one loop we get: 



Ji^ 



(j(i) _ ,(1) _ ^(D) j(o) + ^(1) j(0) ^ [^(1), j(0)] + o(,0) 



(C.2) 



The separation of jet and soft functions is scheme dependent. We work in the formfactor scheme [l3| 
where the jet function is a product of the square root of the formfactors corresponding to each external leg, 
and similarly in the massive case [l5|. It then immediately follows that J*^^^ — i*-^^ = /g /2, where /g ' is 
the one- loop correction to the UV- renormalized gluon form- factor 35|, |36[ : 



/•(I) = ^ 



_^_^ + O(,0) 



(C.3) 



where /3o = llCyi/6 — Np/'i and as is the MS renormalized coupling at scale /x. 
The explicit results for the soft functions read: 



r(l) 



^ i V s T -T 



*7^J = 1 



s: 



(1) 

ah 



= ^^'^'5ab + ^ ^Sg,{Tg)^^^-T,, 



(C.4) 



and the index g denotes the soft-gluon leg. The color matrices pertaining to the soft-gluon are {Tg)cb — ifcab- 
The functions Sij = Sji, not to be confused with partonic invariants, |j depend on whether the legs i,j are 
both massive or not and can be found, for example, in Eq. (29) of Ref. [21| . 

Combining the above results and after some algebra we derive the following expression for the poles of 
the UV renormalized one- loop soft-gluon current: 






1 

2^ 



1 

2^ 






2(pz ■q){pj -q) 



The function hij reads (see, for example, Eq. (29) of Ref. |2l|): 



ln(l 



2x^ 



hi{x) , when rm ^ 0, rrij ^ , 



- h. 



+ 0(6°) . (C.5) 



(C.6) 



and zero otherwise. The variable x is defined in Eq. ((T2]) . The function hij vanishes in the massless limit. 



In deriving the above result we have used color conservation J27=i '^i ~ ^ ^'^'^ ^^^ identity if T^T^ 



-{CaI2)T^, or alternatively, J^*°^ - ^ 



Ca 



Z^i^j=l -^ i -^ j V^i ^jj- 



'In this paper we do not use the notation Sij to denote partonic invariants. 
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Eq. (IC.SP agrees with Ref . [9] . To that end we need to convert the renormahzed couphng to the bare one 
and recah the overall factor of gs in Eq. ([3]) . We also recall the discussion in Section 0] where we explain 
that we work with Uf = ni active flavors and that heavy quar k loops in external gluon fields are exactly 
compensated by the decoupling relation. See also Refs. 137|,[38|, for more information on that point. 

Eq. (IC.5I) applies to space-like kinematics as in Case 1. Continuation to any other kinematics is trivial; 
see Ref. [25| and [Appendix E)| 

Appendix D. Analytical continuation to physical kinematics 

It is often the case that calculations of scattering amplitudes are easier to perform in unphysical kinemat- 
ics. Then the question arises how to analytically continue the result derived in such unphysical kinematics 
to the physical region. In this work, the continuation involves the momentum pj, i.e. we need to continue 
results derived in kinematics where Pj is incoming to kinematics where pj is outgoing. 

The analytical continuation is more involved when pj is massive. When the nromentum pj is incoming 
we have evaluated the soft current exactly in d-dimensions (see Eqns. dH IA.21 IA.41 lA.lll IA.13[ IA.14I) . To 
continue the result to the timelike kinematics where pj is outgoing one has to first express the result in 
a minimal number of variables. Given the scaling-invariance properties of the result, only two variables 
are truly independent. As such we take x and a^, defined in Eqns. (J9I12I) . The variables aj and v can be 
eliminated through the relations aiaj = 1 — t;^ and i; = (1 — a;^)/(l + x'^). 

The rules for the analytical continuation in terms of the variables x and Ui (for general values of the 



masses rriij) are simple, see also Ref. |25| : 



x -^ —X + iS , 

ai ^ at . (D.l) 

The invariance of a^ is easy to understand, since a; ^ {pj ■ q)/{pi ■ Pj)- Not only is a^ invariant under 
Pj -^ —pj □ but, more importantly, its log is: In(ai) — ln((pj • q)) — ln((pi • pj)) -|- . . . = inv. That is 
distinct from the case of aj ^ ^/{{Pj ■ q){Pi ■ Pj)) which itself is invariant but its log is not: In(aj) — 
— ln((pj ■ q)) — ln((pi ■ pj)) + . . . ^ inv. That aj must transform nontrivially also follows from the identity 
aiaj = 1 ~ v^ . 

The practical implementation of the analytical continuation procedure requires the explicit extraction of 
all branching-point singularities around the point a; = 0. Here is a typical example: 

Li2 (-^) = -U2{-x^) - y - 2 \n\x) , (D.2) 

and: 

\n(x^)^\n(x^)+2iTT ; Li„(l - x^) ^ Li„(l ~ x^) - 2i7r^V-^-^^ • (D.3) 

{n-iy. 

We have verified that with the help of Eq. (jD.ll) we can reproduce the first three orders in e of the 
directly calculated integral Mg , from the spacelike calculation of M^ ; see [Appendix A[ Starting from 
the fourth order in e one would have to devise a similar procedure for the set of multiple polylogarithms 
that begin to appear, see Eq. (IA.13|) . This presents an interesting direction for future work that can benefit 
from the number of recent applications of this class of functions in the context of gauge amplitudes in TV = 4 
SYM theories ^. 

When pj is massless, the analytical continuation allows one to obtain the results for configuration Case 2 
from the one for Case 1. For rrij — the conformal variable vanishes, a; = 0, which implies that the analytical 
continuation (|D.1I) becomes trivial. That can also be seen with a direct inspection of the integrals (|A.2IA.3p 
and (|A.7|) in the case rrij — 0. 



^Strictly speaking one inverts not the momentum pj but the sign of all invarians (k ■ pj) linear in pj. 

14 



References 

[1] W. T. Giele, E. W. N. Glover, Phys. Rev. D46, 1980-2010 (1992). 

[2] W. T. Giele, E. W. N. Glover, D. A. Kosower, Nucl. Phys. B403, 633-670 (1993). •hep-ph/ 9302225] . 

[3] S. Frixione, Z. Kunszt, A. Signer, Nucl. Phys. B467, 399-442 (1996). hep-ph/9512328 . 

[4] S. Catani, M. H. Seymour, Nucl. Phys. B485, 291-419 (1997). hep-ph/9605323 . 

[5] S. Frixione, Nucl. Phys. B507, 295-314 (1997). hep-ph/9706545 



[6] S. Catani, S. Dittmaier, M. H. Seymour, Z. Trocsanyi, Nucl. Phys. B627, 189-265 (2002). ' hep-ph/0201036] . 

[7] Z. Bern, V. Del Duca, C. R. Schmidt, Phys. Lett. B445, 168-177 (1998). hep-ph/9810409_. 

[8] Z. Bern, V. Del Duca, W. B. Kilgore, C. R. Schmidt, Phys. Rev. D60, 116001 (1999). |hep-ph/9903516] . 

[9] S. Catani, M. Grazzini, Nucl. Phys. B591, 435- 454 (2000). hep-ph/000 7142^ . 

[10] M. Czakon, Phys. Lett. B693, 259-268 (2010). |arXiv:1005.6 274 [hep-phj]"^ 

[11] M. Czakon, arXiv:1101.0642 [hep-ph]]. 

[12] A. Sen, Phys. Rev. D28, 860 (1983). 

[13] G. F. Sterman, M. E. Tejeda-Yeomans, Phys. Lett. B552, 48-56 (2003). |hep-ph/02 10130]. 

[14] S. M. Aybat, L. J. Dixon and G. Sterman, Phys. Rev. D 74, 074004 (2006) [arXiv:hep-ph/0607309] . 

[15] A. Mitov, S. Moch, JHEP 0705, 001 (2007). hep-ph/0612149 . 

[16] E. Gardi and L. Magnea, JHEP 0903, 079 (2009) arXiv:0901.1091 [hep-ph]]. 

[17] N. Kidonakis, Phys. Rev. Lett. 102, 232003 (2009) arXiv;0903.2561 [hep-ph]]. 

[18] T. Becher and M. Neubert, JHEP 0906, 081 (2009)" arXiv:0903.1126 [hep-ph]]. 

[19] A. Mitov, G. Sterman and L Sung, Phys. Rev. D 79, 094015 (2009) arXiv:0903.3241 [hep-ph]]. 

[20] T. Becher and M. Neubert, Phys. Rev. D 79, 125004 (2009) arXiv:0904.1021 [hcp-ph]]. 

[21] M. Czakon, A. Mitov, G. F. Sterman, Phys. Rev. D80, 074017 (2009). arXiv: 0907. 1790 [hcp-ph]]^ 

[22] A. Ferroglia, M. Neubert, B. D. Pecjak and L. L. Yang, Phys. Rev. Lett. 103, 201601 (2009) ar Xiv:0907.479T1 [hep-phi] . 

[23] A. FerrogUa, M. Neubert, B. D. Pecjak and L. L. Yang, JHEP 0911, 062 (2009) arXiv:0908.3676 [hep-ph]]. 

[24] L. J. Dixon, E. Gardi and L. Magnea, JHEP 1002, 081 (2010) arXiv:0910.3653 [hep-ph]]. 

[25] A. Mitov, G. F. Sterman, L Sung, Phys. Rev. D82, 034020 (2010). arXiv;1005.4646 [hep-ph]]. 

[26] E. Gardi, E. Laenen, G. Stavenga and C. D. White, JHEP 1011, 155 (2010) arXiv:1008.0098 [hep-ph]]. 

[27] A. Mitov, G. Sterman, I. Sung, Phys. Rev. D82, 096010 (2010). arXiv: 1008.0099 [hcp-ph]]. 

[28] J. G. Korner, Z. Merebashvili and M. Rogal, Phys. Rev. D 71, 054028 (2005) .arXiv:hep-ph/0412088] . 

[29] T. Huber and D. Maitre, Comput. Phys. Commun. 175, 122 (2006) arXiv: hep-ph/0507094j ^ 

[30] S. Weinzierl, Comput. Phys. Commun. 145, 357 (2002). math-ph/0201011 . 

[31] S. Moch and P. Uwer, Comput. Phys. Commun. 174, 759 (2006) arXiv:math-ph/0508008 . 

[32] S. Moch, P. Uwer and S. Weinzierl, J. Math. Phys. 43, 3363 (2002) arXiv:hcp-ph/011 0083]. 

[33] J. VoUinga and S. Weinzierl, Comput. Phys. Commun. 167, 177 (2005) arXiv:hep-ph/ 0410259| . 

[34] J. Gluza, A. Mitov, S. Moch et ai, JHEP 0907, 001 (2009). | ,arXiv:0905.1137. [hep-ph]]. 

[35] P. A. Baikov, K. G. Chetyrkin, A. V. Smirnov, V. A. Smirnov and M. Steinhauser, Phys. Rev. Lett. 102, 212002 (2009) 

arXiv:0902.3519 [hep-ph]]. 

[36] T. Gehrmann, E. W. N. Glover, T. Huber, N. Ikizle rli and C. Stude rus, JHEP 1006, 094 (2010) |arXiv: 1004!3653l [hep-phll . 

[37] T. Becher, K. Melnikov, JHEP 0706, 084 (2007). arXiv:0704.3582 [hep-ph]]^ 

[38] M. Czakon, A. Mitov, S. Moch, Nucl. Phys. B798, 210-250 (2008). arXiv:0707^4l39l [hep-ph]] . 

[39] A. B. Goncharov, M. Spradlin, C. Vergu and A. Volovich, Phys. Rev^ Lett. 105, 151605 (2010) ar Xiv:1006.5703l [hep-th11. 



15 



